The nuclear spin-lattice relaxation rate in the PuMGas materials 



Yunkyu Bang 

Department of Physics, Chonnam National University, Kwangju 500-757, 
and Asia Pacific Center for Theoretical Physics, Pohang 790-784, Korea 



O 
O 



o 



c3 



I 

C 

O 

o 



(N 
> 
(N 

O 
m 
o 

O 

03 



i 

C 

O 

o 



X 



M. J. Graf, N. J. Curro, and A. V. Balatsky 

Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
(Dated: February 6, 2008) 

We examine the nuclear spin-lattice relaxation rates 1 /7j of PuRhGas and PuCoGas, both in the supercon- 
ducting and normal states, as well as their Knight shifts. Results for both compounds are consistent with a 
superconducting gap of d-wave symmetry in the presence of strong impurity scattering, though with quite dif- 
ferent gap-over- T c ratios 2Ao/IcbT c (8 for PuCoGas and 5 for PuRhGas). In the normal state, PuRhGas exhibits 
a gradual suppression of (TiT) -1 below 25 K, while measurements for PuCoGas reveal a monotonic increase 
down to T c . We propose that this behavior is consistently understood by the crossover from the two-dimensional 
quantum antiferromagnetic regime of the local 5/-electron spins of Pu to the concomitant formation of the 
fermion pseudogap based on the two-component spin-fermion model. 

PACS numbers: 74.20,74.20-z,74.50 



I. INTRODUCTION 

The discovery of superconductivity (SC) in plutonium 
based systems such as PuCoGas and PuRhGas has stimulated 
the study of unconventional superconductivity and the pair- 
ing symmetry and mechanism in these materials^ The sym- 
metry of an unconventional superconductor is reduced com- 
pared to the symmetry of its normal state, thus resulting in 
many novel properties of the quasiparticle excitation spec- 
trum. With a record high superconducting transition temper- 
ature T c of the order of 20 K among the 4/ and 5/ electron 
based compounds, this class of materials provides hope for 
a unifying pairing mechanism from heavy fermion supercon- 
ductors to the high-T £ cupratesi^i It is believed that the su- 
perconducting action in Pu-1 15 [PuMGas with M=Co and Rh] 
derives itself from the unique character of the 5/ electrons of 
plutonium. 4 PuMGas is isostructural to the tetragonal Ce-1 15 
series [CeMIns]. Very recently, Curro and coworkers^ pro- 
posed, based on their measurements of the Knight shift and 
spin-lattice relaxation rates, that the Pu-1 15 compounds are 
bridging the superconducting and normal-state properties of 
the heavy-fermion Ce- 115 and high-temperature copper-oxide 
superconductors. Therefore providing a means for tuning the 
interaction strength of antiferromagnetic spin fluctuations to 
intermediate values between both extreme limits. 5 

In this paper, we examine the spin-lattice relaxation rates 
Tp 1 of PuCoGas and PuRhGas in the superconducting and 
normal states. Very recently, Sakai and coworkers 6 have mea- 
sured Tf 1 in PuRhGas in the superconducting state and found 
evidence for lines of nodes in the gap function, just as in 
PuCoGas, but with a much reduced superconducting transi- 
tion temperature r c w 8.5 K versus T c w 18.6 K. Both mate- 
rials differ primarily in their lattice parameters, since Rh has 
a slightly larger atomic radius than Co. For example, the ab 
plane and c axis lattice constants are larger in PuRhGas but 
the ratio c/a is larger in PuCoGas, suggesting that PuRhGas 
is less two dimensional (2D). Predicting the consequences 
of this elongation on the superconducting instability is not 



straightforward; however, we find that this small change of 
lattice constants leads to dramatic changes in the electronic 
and magnetic properties in this class of materials. 



Here, we give a detailed theoretical description of the spin- 
lattice relaxation rate and predict what should be observed 
if the Knight shift was measured on the same sample. Our 
self-consistent treatment of impurity scattering in the super- 
conducting state goes beyond the two-fluid approach used by 
Sakai et al. 6 , which was used to explain the large residual 
density of states in PuRhGas. Further, our theoretical fits 
to the experimental data of PuRhGas lead us to three impor- 
tant conclusions: First, the measured pair-breaking effect of 
impurities in PuRhGas reduces the transition temperature of 
a hypothetically pure sample by only 0.5 K. Therefore, the 
lower T c of PuRhGas is an intrinsic property and is not due 
to impurities. Second, the theoretical fits indicate that the ra- 
tio 2Ao /kgT c is markedly reduced in PuRhGas (~ 5) versus 
PuCoGas (~ 8). This fact indicates that the mediating bosonic 
pairing glue is stronger in PuCoGas. Assuming that Pu-1 15 
compounds are spin-fluctuation mediated superconductors, 4 
we conclude that Rh substitution reduces the strength of the 
mediating spin fluctuations. Indeed, this conclusion is sup- 
ported by the observed behavior of 1/7] in the normal state 
of both compounds. Third and most interestingly, when the 
experimental data are plotted as T\ T versus T, we find that in 
the normal state T\ T of PuRhGas saturates and is nearly flat 
over a wide temperature region T c < T < 3T C , whereas T\T 
in PuCoGas shows a monotonic decrease down to T c . This 
saturation resembles closely the pseudogap (PG) feature of 
the underdoped cuprates. 7 As a possible explanation of this 
phenomenon, we propose the two-component spin-fermion 
model of antiferromagnetically correlated metals 8 and argue 
that the two natural energy gaps accounting for the PG behav- 
ior in Pu-1 15 are the spin gap A$g and the fermion gap Ape. 
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II. SUPERCONDUCTING STATE 

The experimental techniques of nuclear magnetic reso- 
nance (NMR) and nuclear quadrupolar resonance (NQR) have 
been used successfully in the past to distinguish between the 
spin states of Cooper pairs (spin singlet vs. spin triplet pair- 
ing) and provide indirect information on the symmetry of the 
gap function - fully gapped vs. nodal lines or nodal points in 
the gap function on the Fermi surface. Both techniques probe 
directly the quasiparticle density of states and reveal indirect 
information about the pairing symmetry. 

The standard explanation of power vs. exponential laws in 
the low-temperature behavior of thermodynamic and transport 
properties, for example, the spin-lattice relaxation rate l/7\, 
comes from the difference of nodal and fully gapped excita- 
tion spectra in the superconducting state. In clean nodal su- 
perconductors l/T\ exhibits a nearly T 3 behavior far below 
the superconducting transition temperature, T <C T c , while it 
is exponential for gapped superconductors. On the other side, 
deviations from this behavior, like the ^-linear temperature 
dependence of 1 /T\ at low temperatures, are explained by im- 
purity effects in an unconventional superconductor with lines 
of nodes on the Fermi surface. 

In our calculations the effect of impurity scattering is in- 
cluded within the self-consistent T-matrix approximation, 9 
which is the standard formulation for pointlike defects in a su- 
perconducting dilute alloy JSiiiii^ For the case of particle-hole 
symmetry of the quasiparticle excitation spectrum the Nambu 
component 73 of the T matrix vanishes, and for a d-wave order 
parameter with isotropic scattering T\ = (also without loss 
of generality we can choose T% = by general U(l) gauge 
symmetry), where 7; is the rth component of the 2x2 Nambu 
matrix expanded in Pauli matrices. Then we need to calculate 
only 7b (co). The impurity self-energy is given by Eo = FTq, 
where F = nj/nNo- Here No is the normal density of states 
(DOS) at the Fermi surface, is the impurity concentration; 
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T °«) = [J-iL)] > where So(««) = m- ^k a,2 +£ i +A 2 W - 
impurity renormalized Matsubara frequency is defined by 
CO,, = co„ + Eo, with co„ = nT(2n + 1), and the scattering 
strength parameter c is related to the s-wave phase shift 8o 
by c = cot(8o). Using this self-energy Eo the following gap 
equation is solved self-consistently, 
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where V(cj> — ((>') is the angular parameterization of the pairing 
interaction, and co c is a typical cutoff energy. We assume the 
canonical d-wave gap function of the form A(k) = Ao(cos k x — 
cos ky) or A((|)) = Aocos(2(j)) for a cylindrical Fermi surface. 
The pairing potential V ((]) — (])' ) induces a gap with d-wave 
symmetry. Although its microscopic origin is not the issue of 
this paper, we believe it originates from 2D antiferromagnetic 
(AFM) spin fluctuations. The static limit of the spin suscep- 



tibility of the AFM fluctuations, %(q, CO = 0) 
parameterized near the AFM wave vector Q as 8 



V($-$)=V d (b)- 
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where the parameter b is inverse proportional to the AFM cor- 
relation length ^, normalized by the cylindrical Fermi surface 
(£, ~ an/b; a is the lattice parameter). For all calculations in 
this paper, we choose b = 0.5 which is not a sensitive parame- 
ter for our results unless ^ is very large (b < 0.1), 8 i.e., within 
the range of 0.1 < b < 1 our results show little variations and 
are qualitatively the same. 

With the gap function A((j)) and Tq (co) obtained from Eq. Q 
(where Tq (co) is analytically continued from 7b(co„) by Pade 
approximation) we calculate the nuclear spin-lattice relax- 
ation rate 1/Ti following the standard formulation^2ii£ 
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(3) 

and the superconducting spin susceptibility 
Is f°°df FD ((£>) /_ co 




3co 



Re 



(4) 



where fp d(co) is the Fermi-Dirac function, the impurity renor- 
malized quasiparticle energy CO = co + Eo(co), and (...)§ means 
the angular average over the Fermi surface. The first term in 
the bracket of Eq. is N 2 ((ii). The second term vanishes 
in our calculations because of the symmetry of the gap func- 
tion. To calculate \/T\T using Eq. or Xs using Eq. @, 
we need the full temperature dependent gap function A(0, T) 
and T c , Our gap equation Eq. Q is the BCS gap equation, 
therefore it gives the BCS temperature behavior for A(<j),r) 
and Ao = 2.14£g71- for the standard weak-coupling d-wave 
SC. In order to account for strong-coupling effects we use 
the phenomenological formula A(ty,T) = A(§,T = 0) Z(T) 
with E ( T ) = tanh ( P ^/T c /T — 1) , and parameters p and Ao / T c . 
Then we only need to calculate A(0,O) at zero tempera- 
ture. The temperature dependence of Eo(co,r) = r7b(co,r) 
is similarly extrapolated: 7b(co, T) = Tq(q),T = 0) E(T) + 
T nor mal(l ~ E(T)), where T norma i = F/(c 2 + 1) is the normal 
state 7b. In our numerical calculations we chose p = 1.74, 
because our final results are not very sensitive with respect to 
this parameter, while the ratio Ao / kgT c is an important param- 
eter to simulate strong-coupling effects. The larger the gap 
ratio the more important are strong-coupling effects. 

In figures [2 and |2] the spin-lattice relaxation rate of 
PuRhGas by Sakai et al. 6 is shown. For ease of comparison 
l/7i was normalized to l/7i = 10 at T = T t .. The insets show 
the corresponding normalized quasiparticle DOS for varying 
scattering rates F. With our earlier described choice of gap pa- 
rameters (p = 1.74 and 2Aq = 5ksT c ) an impurity scattering 
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FIG. 1: The NQR spin-lattice relaxation rate 6 plotted versus tem- 
perature normalized by the superconducting transition temperature 
T c = 7.6 K. Calculations are shown for 2Aq = 5ksT c and three val- 
ues of the impurity scattering rate T for unitary scattering. Inset: The 
normalized quasiparticle density of states is shown for corresponding 
values of the impurity scattering rate T/Aq = 0,0.032,0.064. 




FIG. 2: The NQR spin-lattice relaxation rate 6 plotted versus tem- 
perature normalized by the superconducting transition temperature 
T c = 7.6 K. Calculations are shown for 2Ao = SksT c and three val- 
ues of the impurity scattering rate T for unitary scattering. Inset: The 
normalized quasiparticle density of states is shown for corresponding 
values of the impurity scattering rate T/Aq = 0,0.032,0.064. 



rate of T/Aq = 0.032 in the unitary limit (c = 0) is enough to 
completely fill the low energy gap with impurity states, where 
N((0 = 0) reaches more than 25% of the normal-state DOS 
No. For a higher value of T/Ao = 0.064, the T-linear region 
extends up to ~ 0.35 T c . At temperatures near T c the coher- 
ence peak is almost invisible because of the sign-changing gap 



FIG. 3: The calculated spin susceptibility Xs °f a d-wave SC nor- 
malized by its normal state value %N f° r g a P values 2Ao = 5kgT c 
(solid lines) and 8£#r c (dotted lines), and impurity scattering rates 
r/Ao = and 0.032. 

function, i.e., vanishing of the second term in Eq. 0. 

Below T c the spin-lattice relaxation rate shows a nearly T 3 
behavior until it crosses over to a T-linear region. A com- 
parison with the experimental data by Sakai et al. 6 is in good 
agreement with a scattering rate close to T/Aq = 0.032. Based 
on this value, we estimate the reduction of T c q of a hypotheti- 
cally pure sample due to impurities to be T c o — T c = |T ss 0.53 
K, which results in r c o « 8.1 K or 9.0 K, depending on the 
value of T c = 7.6 K or 8.5 K. This small suppression and large 
impurity induced DOS is a characteristic of unitary impurities 
(c=0). 

In Fig.^ we obtained a better fit to the experimental data 
(symbols) assuming a slightly lower superconducting transi- 
tion temperature T c = 7.6 K than the reported value of T c = 
8.5 K by Sakai et al. 6 This could indicate the presence of 
a pseudogap similar to the high-temperature superconductor 
YBa2Cu307_5, where 1 /T\ is suppressed starting above T c . 

In Fig. |3 we plot 1 /T\ for an enhanced strong-coupling 
d-wave gap of 2Ao = 8fc#r c , as was recently reported for 
PuCoGag. 3 Due to the larger gap value, the theoretical l/7\ 
drops initially faster below T c than the experimental data. 
Hence, we find a poorer fit to the measured data for this choice 
of strong-coupling gap. It demonstrates that the comparison 
of 1 /T\ data with theoretical calculations is a useful tool for 
determining the strong-coupling gap value Ao/£gr c . 

Fig- E] shows the prediction for the spin susceptibility, 5(5, 
or its corresponding NMR Knight shift, K = Kq + A^s, where 
Kq and A are constants for most materials. Xs is calculated for 
the same d-wave gap values as was used for the spin-lattice 
relaxation rates in figures ^ an d|2] Again a modest impurity 
scattering rate of T/Aq = 0.032 results in a large residual sus- 
ceptibility at zero temperature, equivalent to roughly 25% of 
the normal-state DOS or spin susceptibility %n. The quantita- 
tive difference in the spin susceptibilities between gap values 
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2Ao = 5kgT c and 8ksT c should be easily discernible in future 
measurements of the Knight shift. 

At first sight it might appear that the sample PuRhGas had 
three times more defects than PuCoGas of similar age, based 
on our best-fit scattering rates of T/Ao = 0.032 for PuRhGas 
and r/Ao = 0.01 for PuCoGas. This could be explained by 
slightly different isotope mixes of plutonium used. How- 
ever, if we express the scattering rates in absolute values, we 
find very similar impurity scattering rates for both samples, 
namely, r = 0.6 K for PuRhGa 5 and T = 0.7 K for PuCoGa 5 , 
consistent with a common origin of defect generation due to 
self-irradiation by plutonium. 

One final remark is warranted, namely, that the experi- 
mental data by Sakai et al. in Fig.Q are normalized assum- 
ing a slightly lower superconducting transition temperature 
r t = 7.6 K than the reported value by the authors. Indeed, 
if T was normalized by T c = 8.5 K, then the fitting would be 
poorer and in particular the excess relaxation rate just below 
T c could not be explained by a simple superconducting tran- 
sition. The origin of this ambiguity of T c is not clear at the 
moment. If the true T c of the sample was indeed 7.6 K, as was 
used in Fig.^ then the incorrectly assigned T c = 8.5 K might 
be due to the presence of a pseudogap that will be discussed 
in the next section. 



III. NORMAL STATE 

Next we address the non-Fermi liquid behavior of the Pu- 
115 compounds in the normal state. As we will argue be- 
low, this behavior can be consistently understood within spin 
fluctuation theory. In Fig.|4ja), we plot the inverse of the mea- 
sured 1 / T\ for PuCoGas and PuRhGas multiplied by T versus 
temperature up to higher temperatures (multiple times T c ). A 
Fermi liquid should exhibit a constant T\ T in the normal state, 
in contrast to the succinct features of Pu-1 15: (1) at high tem- 
peratures (T > 25 K) both data show T-linear behavior; (2) 
T\T in PuCoGas shows monotonic decrease down to T c — 
a small deviation from the T-linear behavior in the region of 
T c ^ T < 30 K calls for additional explanation; (3) the most 
interesting feature of the curve is the gradual round-off of T\ T 
in the PuRhGas data below ~ 20 K. Even the superconduct- 
ing transition at T c = 8.5 K or 7.6 K is not clearly discernible 
from these data. This roundoff in \/T\T starting far above T c 
has been frequently observed in underdoped high temperature 
superconductors and has been attributed to the suppression 
of low-energy spin fluctuations associated with the pseudo- 
gap behavior. Recently, experiments of several heavy fermion 
compounds showed such a pseudogap behaviorpi^ii While the 
pseudogap behavior in heavy fermions typically shows only in 
a very narrow temperature range of a few Kelvin, it provides 
a much clearer evidence for its magnetic origin than in the 
cuprates. For the cases of Celn3 and CeRhln^i^i^ the pseu- 
dogap occurs in NMR data above the Neel temperature T^. 
In addition, the pseudogap of CeCoIns occurs above the su- 
perconducting transition, a system which is known to be very 
close to the two-dimensional antiferromagnetic quantum crit- 
icality (QC)iii 



Considering that Pu-1 15 materials are near a 2D antifer- 
romagnetic (AFM) instability, we start with the phenomeno- 
logical model of the antiferromagnetic ally correlated metal. 
The minimal set of low energy degrees of freedom are the 
fermionic charge excitations and the collective spin excita- 
tions. This phenomenological theory is also called spin- 
fermion model and has been intensively studied by Pines and 
coworkers^ In contrast to this standard spin-fermion model, 
we proposed for heavy fermions the two-component spin- 
fermion model, 8 where the spin modes originate directly from 
localized spins rather than from collective particle-hole exci- 
tations. In a mixed momentum and real-space representation 
the corresponding Hamiltonian is written as 

H = £4(k)e(k)c a (k) + I /S(r)-4(r)a a pc p (r)+# s 

k,a r,a,p 

(5) 

where the first term is the fermionic kinetic energy and the 
second term describes the coupling between local spins S(r) 
and the spin density of the conduction electrons. The last term 
represents an effective low-energy Hamiltonian for the local 
spins. When the local spins have a short range AFM correla- 
tion, the spin correlation function has the general form 17 



x(q,») = 



x(Q,o) 



1 + ^ 2 |q - Q| 2 - co 2 /A 2 G - /co/co ' 



(6) 



where A$c is the spin gap, Q the 2D AFM ordering vec- 
tor, \ the magnetic correlation length, and CO the spin relax- 
ation energy scale, which comes from Landau damping of 
the fermionic sector. Given the above form of the spin sus- 
ceptibility %(q,co) and assuming the 2D AFM correlation 18 , 
it has been shown that l/TiT.~ 1/cOpi^ Further, it is known 
that CO ~ for the z=l quantum-critical phase of the 2D 
quantum antiferromagnet. 20 As was shown by Chakravarty 
and coworkers, 20 the magnetic correlation length displays, de- 
pending on the phases of quantum criticality (QC) or quantum 
disorder (QD) the following behavior: 



:-\ 



T forQC (T>T*) 
const, for QD (T < T*) 



(7) 



The high-temperature T-linear behavior in T\T, see 
Fig. |3 a), of PuCoGas and PuRhGas can be understood as 
a generic feature of the quantum critical regime of the 2D 
Heisenberg antiferromagneti^P. In principle, this high tem- 
perature QC explanation can be tested by the measurement 
of the spin-spin relaxation time Tig, as was done for the 
high-temperature superconductors; 21 then the simple relation 
T\T IT^q =const. should characterize the QC dynamics in Pu- 
1 15 by determining the dynamic critical exponent z. Also the 
low-temperature round-off in T\ T of PuRhGas can be associ- 
ated with the QC to QD crossover of the 2D AFM. However, it 
is still not sufficient to explain the additional fall-off of \/T{T 
between 20 K and T c w 7.6 K (for more detail see the inset of 
Fig. 0b)). We propose that this additional suppression of the 
spin-fluctuations is caused by the suppression of the fermionic 
DOS (namely the fermionic pseudogap). Since the term /co/co 
in Eq.(6) originates from Landau damping of the fermionic 
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FIG. 4: (a) Plot of T\T data vs. T of PuCoGag and PuRhGa 5 . (b) 
Plot of l/T\ T vs. T for a wider temperature range of the normal state 
for PuRhGa5 (open diamonds) and PuCoGas (open circles) and the 
theoretical fits. Parameters for the theoretical fits are A^c =20 K and 
AppQ =8 K for PuRhGas and A$c =0 K and Appc =8 K, respectively. 
Inset: Close-up view for PuRhGa5 at low temperatures. Theoretical 
fits with (dotted line) and without (solid line) the fermionic PG cor- 
rection are shown for comparison. (PuRhGas data are from Ref.[6].) 



sector, instead one needs to include ~ mN(Ep ,T) /a>. There- 
fore, when the fermionic DOS N(Ep) becomes temperature 
dependent, it needs to be included and leads to the modifica- 
tion 1/TiT ~N(E F ,T)/G)(T). 

Several authors^ have studied the influence of magnetic 
correlations on fermionic quasiparticles and found that in- 
creasing the magnetic correlation length Z, causes a precursor 
effect of a spin-density wave state, which forms a quasi-gap in 
the DOS. The amount of the suppression of the DOS depends 
sensitively on the parameters, such as the size of the corre- 
lation length the coupling constant, and temperature, etc. 



In this paper, therefore, we merely introduce a phenomeno- 
logical form of the pseudo-gapped DOS, N(Ef,T,Afpg). 
The systematic numerical studies of the fermionic pseudo- 
gap (FPG) Afpg due to magnetic correlation in the two- 
component spin-fermion model will be reported elsewhere. 

Combining the temperature dependence of the FPG and the 
magnetic correlation of the 2D AFM, we can write \/T\T in 
the two-component spin-fermion model as 



\/T{T = AjV (£/0[l-tanh 2 ( 



&FPG 



2VT 2 +r 2 ' 



B. 



(8) 



[A 5G + rexp(-4A SG /r)] 



The first factor No(Ep)[l ~ ■■■} is the phenomenological 
form of the fermionic DOS with FPG Appc and damping rate 
r. The second factor tt— — , is a smooth crossover function 20 

[Asg+-] 

for Z,(T) describing the QC to QD behavior of Eq.(7) with the 
spin gap Asg ~ T* . Here, B is a constant describing a tem- 
perature independent contribution and A is a constant scale 
factor. Using this formula, we fit the experimental data of 
PuCoGag and PuRhGas for normal state in Fig.@Ib). The fit- 
ted results are in excellent agreement with experiment, where 
the two key fitting parameters provide estimates for the im- 
portant energy scales of this phenomenological model. For 
PuRhGas, we used Asg — 20 K and Appc = 8K; the damp- 
ing rate F is not a very sensitive model parameter, so we use 
r=25 K in all cases. A spin gap of Asg ~ T* = 20 K can be 
read off from the data in Fig. [3 a), where T\T deviates from 
a linear temperature dependence. However, notice that with- 
out the FPG correction Afpg the additional drop of \/T\T 
in the region of T c < T < T* cannot be explained only by 
the QC to QD crossover of cj, see the inset of Fig.|4]h). For 
PuCoGas, the monotonically increasing \/T\T at lower tem- 
peratures implies increasing t,(T) and stronger magnetic cor- 
relations than in PuRhGas. Along the same line of thought, 
as applied to PuRhGas, the FPG should also be formed in 
the case of PuCoGas. In Fig.|4jb), we used Asg = K and 
&fpg = 8 K for PuCoGas. The FPG effect is not much visible 
because it is overwhelmed by the stronger temperature depen- 
dence of 1 /co j for PuCoGas. For PuCoGas, we could ob- 
tain a same quality of good fits with different values of Appc 
as large as ~ 20 K. 

We summarize the effects of the magnetic correlations in 
PuCoGas and PuRhGas on the thermodynamic behavior in the 
schematic phase diagram in Fig. [5] with respect to a generic 
coupling parameter g. The assignment of PuCoGas rather 
than PuRhGas to the vicinity of the quantum critical point 
(QCP) in the QC region is consistent with various other ex- 
periments: for instance, the lower T c , the smaller the value of 
the specific heat jump AC/T C w 45 mJ mol -1 K" 2 ^ and the 
smaller the gap-over- T c ratio 2Aq/T c in PuRhGas compared to 
PuCoGas. Also the resistivity of PuCoGas shows an anoma- 
lous power law dependence 1 , p °< T 4 / 3 , as expected near a 
quantum phase transition. All these experimental data indi- 
cate that the 2D AFM fluctuations are weaker in PuRhGas 
than in PuCoGas. In Fig.|5]we added the crossover line (thin 
dotted line) to the FPG region, Tfpg- The formation of the 
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FIG. 5: The schematic phase diagram for the alloy system 
Pu(Co,Rh)Ga5. Thick solid lines are the 2D AFM transition tem- 
perature (TJv) and the QC to QD crossover temperature (T*). The 
thin dotted line is the crossover temperature (Tfpc) to the fermionic 
pseudogap (FPG) region. 



FPG is not a universal property and depends on material spe- 
cific details. Nevertheless, in order to understand the detailed 
behaviors of \JT\T in PuMGas materials, it is indispensable 
to be included. 

To complete our discussion, we briefly describe possibil- 
ity of alternative explanations. In recent Knight shift (K) 
measurements in PuRhGa5, Sakai et ali^i, found that the 
Knight shift failed to track the bulk susceptibility (%) below 
a temperature 30K. Such K-% anomaly is known to exist 
in other heavy fermion compounds and various explanations 
were proposed 25 . Naturally this anomaly is likely to be related 
to the anomalous temperature dependence of Ti— ; however, 
the deviation of K=< % relation perse doesn't necessarily mean 
the suppression of \/T\T . There are three proposals for K-% 
anomaly in the literature: (1) CEF (crystal effective fields) 27 : 
the population of different CEF levels of f-electron changes 
with temperature and hence the HF (hyperfine) coupling be- 
tween nuclei spin and the f-electron spin obtains a tempera- 
ture dependence. (2) Kondo cloud screening 2 -: the HF cou- 
pling between nuclei spin and the f-electron spin is mediated 
by the RKKY modulation of the conduction electrons. The 
onset of Kondo screening would change the characteristics of 
the RKKY modulation and leads to the change of HF cou- 
pling. (3) Two fluids model 25 : total susceptibility consists of 
two part - one from f-electrons and the other from conduction 
electrons - and each component has different temperature de- 
pendence and HF couplings to the nuclei spin. Due to the 
different temperature dependence of each susceptibility, the 
non-proportionality between K and % can be explained. 

The proposals (1) and (2) are basically invoking on the tem- 
perature dependent HF coupling but by different mechanisms: 
Kondo and CEF, respectively. Our model is similar to the pro- 



posal (3) in the spirit of the two fluids model. However, there 
are technical and conceptual differences between (3) and our 
model. First, technical difference is the following. While we 
also assume two susceptibilities - one from local f-electron 
and the other from conduction electron, our two fluids model 
is not a simple addition of two susceptibilities. We assume 
that the dominant contribution ( the interesting temperature 
dependent part) always comes from the local f-electrons. But 
the f-electron susceptibility essentially obtains its low energy 
dynamics through coupling with the conduction electrons and 
as a result the conduction electron susceptibility always feeds 
back into the f-electron susceptibility and vice versa (see /co/co 
term in Eq.(6)). Second, the most important conceptual dif- 
ference is that we assume that the temperature dependence of 
the f-electron susceptibility for the range of interesting tem- 
peratures is arising from a magnetic correlation of local f- 
electrons but not from Kondo or CEF effect. This point of 
departure is very crucial and should be determined by exper- 
iments. There are already abundant experimental evidences 
that CeM (M =Co, Rh, Ir)Ins materials are inside or in the 
vicinity of AFM ordered phase 29 . For PuM (M =Co, Rh)Ga 5 , 
there are yet no direct measurements of the magnetic corre- 
lations, but the fact of a D-wave pairing 3 - in the supercon- 
ducting state and the anomalous temperature dependence of 
dc-resistivity 4 indicate that PuMGas materials are near AFM 
instability as sketched in Fig. 5. 



IV. CONCLUSION 

We have studied the nuclear spin-lattice relaxation rates 
l/7i of PuCoGas and PuRhGas, both in the normal and su- 
perconducting states. In the superconducting state, both com- 
pounds display the features of a dirty d-wave superconduc- 
tor with impurity scattering in the unitary limit. This also is 
borne out in the calculated Knight shifts. The superconduct- 
ing gap values of 2Ao/ T c are ~ 8 and ~ 5 for PuCoGas and 
PuRhGas, respectively, indicating that pairing fluctuations are 
much stronger in PuCoGas than in PuRhGas. 

In the normal state, the temperature behavior of 1 / T\ be- 
tween both compounds is qualitatively different at low tem- 
peratures. While \/T\T of PuCoGas displays a genuine 
quantum-critical behavior similar to a 2D quantum antiferro- 
magnet down to T c , PuRhGas shows a pseudogap-like sup- 
pression over a wide temperature region from T c (7.6 K) to 
roughly 3T C (25 K). Because of this remarkable observation, 
we proposed the two-component spin-fermion model and ar- 
gued that the magnetic spin gap originates on the local spins of 
the 5/ electrons of Pu and that the concomitant formation of 
the FPG (Afpo) can provide a consistent explanation of these 
phenomena. 

We argued previously 3 - 4 that PuCoGas bridges the heavy 
fermion superconductors and high-r c cuprates in terms of 
the superconducting pairing mechanism. The observation of 
the pseudogap phenomenon in PuRhGas and its temperature 
range of over two times T c (or roughly 15 K) is another evi- 
dence that PuMGas is indeed the missing link between heavy 
fermion superconductors and cuprates, holding the key in- 
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gredient - magnetic correlations - yet with an intermediate 
energy scale. To test and confirm this hypothesis, we pro- 
pose systematic studies of the alloy system Pu(Co,Rh)Gas. It 
will be an ideal system for further studies because it can be 
cleanly tuned to explore the magnetic phase diagram sketched 
in Fig.|5]without changing the carrier density in contrast to the 
high-r c cuprates. 
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